Instanton Theory of Burgers Shocks and Intermittency 



OO 
O 

o 

O 

Q 

OS 



c : 

■ 

i < 
3 ■ 

o ■ 

C/J ■ 

>v 
^ : 

Oh. 



(N 
00 

o 

> 

X 
S3 



L. Moriconi 

Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
CP. 68528, 21945-970, Rio de Janeiro, RJ, Brazil 

A lagrangian approach to Burgers turbulence is carried out along the lines of the field theoretical 
Martin-Siggia-Rose formalism of stochastic hydrodynamics. We derive, from an analysis based on 
the hypothesis of unbroken galilean invariance, the asymptotic form of the probability distribu- 
tion function of negative velocity-differences. The origin of Burgers intermittency is found to rely 
on the dynamical coupling between shocks, identified to instantons, and non-coherent background 
fluctuations, which, then, cannot be discarded in a consistent statistical description of the flow. 
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I. INTRODUCTION 

The long established Burgers model of compressible 
one-dimensional flow [l[ provides an interesting testing 
ground for the performance of numerical and analytical 
strategies in turbulence research. Despite its own pecu- 
liar phenomenology, as evidenced by the complete fail- 
ure of approximations based on the "K41" scaling theory 
, there are important conceptual analogies between the 
Burgers model and usual three-dimensional turbulence 
(see Ref. Q for a comprehensive review). 

We note, in passing, that the Burgers model is more 
than just a "mathematical toy"; in its multidimensional 
version the Burgers equation plays an important role in 
several realistic problems, like nonlinear acoustics 
cosmology [3, critical interface growth [?J , and traffic 
flow dynamics [q. 

A great deal of attention has been focused on the prob- 
lem of non-gaussian fluctuations observed in the high 
Reynolds number regime of Burgers dynamics - the inter- 
mittency phenomenon, for short. As it is verified through 
numerical simulations [ic| , velocity differences 



z = u(x + C, t) — u{x — C, t) 



(1.1) 



are found to be very intermittent at small scales (C much 
lesser than the integral scale L). The probability distri- 
bution function (pdf) of the right tail of z, which decays 
faster than gaussian, has been analytically obtained in a 
number of different ways 0, HO, El ■ The left tail, 
on the other hand, which is related to Burgers shocks, is 
found to have a power law profile p(z) ~ l/|z| Q , with no 
sharp consensus on the value taken by the exponent a. 

A Fokker-Planck approach to the computation of 
velocity-difference pdfs, with closure given by an operator 
product expansion treatment of the dissipative anomaly 
was put forward by Polyakov [l2j ■ This method provides 
a fine description of the pdfs right tail, and yields a 
power law form for the left tail with 5/2 < a < 3 fl4| . 
Extensive numerical simulations performed by Gotoh and 
Kraichnan (Io| indicate that a — 3. At variance with 
such findings, an analytical study based on the velocity 
field profiles in space-time neighborhoods of shocks, the 
so-called preshock events, gives a = 7/2 (TJ, fl5f . a re- 



sult confirmed by alternative lagrangian simulations of 
the Burgers equation [l?], EH • 

In an attempt to conceal these apparently contradic- 
tory conclusions, Boldyrev et al. have suggested that the 
left tail exponent is not universal, departing from a = 3 if 
flow realizations fail to satisfy a strong form of galilean in- 
variance fl9( | . which holds, by definition, if usual galilean 
invariance is observed in the bulk, regardless the bound- 
ary conditions at infinity. In rephrased form, the whole 
point of Ref. [l9[ is that finite-size effects which break 
strong galilean invariance would lock larger fluctuations 
of shock jumps and negative velocity derivatives, reduc- 
ing intermittency. In this paper, we find support to the 
conjecture that the left tail exponent is a = 3 when the 
strong form of galilean invariance is fulfilled. 

This work is organized as follows. In section II, we in- 
troduce Burgers intermittency as a phenomenon related, 
in the inviscid limit, to shock amplitude fluctuations. 
The great convenience of a lagrangian description of the 
flow is then pointed out. In section III, we discuss the 
Martin-Siggia-RosefMSR) formulation of stochastic hy- 
drodynamics [2(| [2l| within the lagrangian perspective. 
In section IV, Burgers shocks will be given as instantons 
[Til fl6| , and background fluctuations around them will 
be taken into account in the computation of the asymp- 
totic behaviour of the pdf of negative velocity differences. 
In section V, we summarize and discuss our results. 

In appendices A and B we provide technical details on 
some of the material discussed in section IV. 



II. PHENOMENOLOGICAL CONSIDERATIONS 

The Burgers model describes the dynamics of a one- 
dimensional velocity field, u = u(x,t), ruled by the evo- 
lution equation 



dtU + ud x u — vd\u + / 



(2.1) 



where v is the kinematical viscosity and / = fix, t) is the 
external force which sustains the flow and introduces the 
integral length scale L. There is no pressure term in the 
above equation, and no imposition of incompressibility 
as well. 
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Let uq(x) = u(x, 0) be the velocity field at initial time, 
supposed to be a C 1 function defined on — oo < x < oo. 
The Cauchy problem is exactly solvable for the Burgers 
equation [22, l23|. In the forceless case, the velocity field 
is, at time t > 0, 



u(x,t) = —2vd x \n 



dy exp 



1 f v 

— / dxuo(x ) 

£V In 



(2.2) 



As it is well-known, (|2.2|) leads, in the vanishing viscosity 
limit, to discontinuous shocks (i.e., the velocity field be- 
comes piecewise C 1 ), which can be interpreted as "sinks" 
of fluid particles. External forcing does not spoil the 
process of shock creation, even though it can affect the 
statistics of shock amplitudes. 

Stable and unstable regions of the flow are distin- 
guished, essentially, by the sign of the spatial velocity 
derivative. Neglecting higher order corrections, consider 
the expansion u(x,t) — cr (<) + a\{t)x in the neigh- 
borhood of an arbitrary point. Eq. (|2.1j) gives, for 
f(x,t)=0, 



(T + CTQCTi = 

oi + o\ = , 



leading to 



where 



cro(t) = <ro(0) exp 



o-i (t) = 



dt'ai{t') 



«7!(0) 



l + o-i(0)t 



(2.3) 



(2.4) 



(2.5) 



Therefore, if oi(0) = d x u\t=o is positive, we expect, from 
(|2.5j) . that d x u will decay. On the other hand, if o-i(O) is 
negative, then \d x u\ will increase in time, implying flow 
instability. This is the mechanism for the generation of 
large negative velocity derivatives in Burgers turbulence, 
which in a time scale of order l/|9 x w| are transformed 
into long lived shocks. 

We are interested to study the statistics of negative 
velocity differences z for v — > and Q/L <C 1. Under 
these circumstances, shock jumps provide the main con- 
tribution to the strong fluctuations of z. A time series of 
z would exhibit intermittent negative spikes associated 
with the passage of shocks, separated by much weaker 
signals due to smooth velocity configurations. 

Suppose, now, that fluctuations of z are alternatively 
measured from the subtraction of velocity fields defined 
at points x(t) + C and x(t) — £, where x(t) is the po- 
sition of a fluid element that moves with the flow. We 
expect to have, in such a lagrangian framework, the same 
asymptotic power law form for the left tail pdf of z. The 
central point in this correspondence is that fluid particles 
typically spend a finite fraction of their times at shock 



discontinuities. Once a fluid particle is dragged into a 
shock discontinuity, it remains there until the shock col- 
lapses, or it is absorbed by another one. 

Negative spikes in the eulerian time series of z are re- 
placed, in the lagrangian reference frame, by smooth fluc- 
tuations of shock amplitudes, which last for much longer 
times. In order to compute statistical properties of the 
eulerian negative spikes, one would have to find out how 
shocks match to each other in solutions of the Burgers 
equation. Within the lagrangian framework, on the other 
hand, it suffices to describe fluctuations around isolated 
shocks - a much simpler task that points out, for our pur- 
poses, the advantage of lagrangian methods over eulerian 
ones. 



III. PATH-INTEGRAL FRAMEWORK 

In the stochastic hydrodynamics approach to Burgers 
turbulence, the external forcing in (|2.1[) is taken to be 
a large scale gaussian random field, with zero mean and 
correlator 

(f(x,t)f(x',t'))=D(\x-x'\)8(t-t') , (3.1) 

where we take 

D(\x-x'\) = D aip(-\x-x'\ 2 /L 2 ) . (3.2) 

The conditional probability density functional to have 
velocity configuration uq(x) at time t = if U-t{x) is 
the velocity at time t = — T, can be written as the path 
integral [Ifl, [H[ 



Z = Af J DuDuexp(iS) 



(3.3) 



where Af is an unimportant normalization factor (which 
will be usually supressed from notation) and 

S = S[u, u] = J dt J dxu[dtu + ud x u — vd x u] 
+ ~2 f dt J dxdx ' {l ( x ' t ) {l ( x '' t ) D (\ x ~ x '\) (3- 4 ) 



is the so-called MSR action. Expressions (|3 - 3(1 and (|3.4|) 
are subject to the boundary conditions 



uo(x) = u(x,0) 
u-t{x) = u{x,—T) 



(3.5) 



In order to pave the way for a lagrangian description 
of the flow, let us consider a general reference frame R' 
which moves with velocity 4>{t) relative to the original 
(inertial) "laboratory" frame R. The position and veloc- 
ity in R' are 



x' = x — / dt'(f>(t') , 



-T 



Uj,(x',t) = u(x,t) - (f>{t) 



(3.6) 
(3.7) 



3 



In the moving frame, the velocity at the origin (x 1 = 0) 
is defined as 



U(f> {t) ee u (O, t ) = u [J *w)> * ) " <K*) ( 3 - 8 ) 

For a given field m = u(x, t), there is a unique time- 
dependent function (f>(t) which solves u^,(0,t) — 0. It 
is clear that <j>(t) is in this case the velocity of a locally 
comoving reference frame - that's how lagrangian coordi- 
nates come into play. We introduce, correspondingly, the 
Faddeev-Popov determinant [24J, A[u(0, t)], by means of 



A _1 [w(0,t)] ee / D#[u (t)] 



(3.9) 



Note that A[u(0,i)] is invariant under the generalized 
galilean tranformations given by (|3.6|) and (|3.7|) . In fact, 



a- 1 [«*,(*)] = y ^^k+0oW] 

= J D(j)5[u^t)] = A- x K0,t)] . (3.10) 
Relation (13. 9p yields 

A[u(0,t)] / WM*)] = 1 . (3.11) 



Inserting (I3.11|) into the integrand of (13. 3|) and exchang- 
ing the order of integrations, we get 



Dcj> I DuDuA[u{0,t)}S[u^{t)}exp{iS) . (3.12) 



Generalized galilean transformations can be used 
to replace the Dirac's delta functional in (|3.12[) by 
<5[u0 =o (i)] = 6[u(0,i)]. To accomplish that, we first sub- 
stitute, in the MSR action of (|3.12j) . the integration fields 
u(x,t) and u{x,t) by galilean transformed ones, through 

u(x, t) — u<j,(x', t) + cj)(t) , (3.13) 
u(x,t) =u ( p(x',t) . (3-14) 



We find 



where 



S = Scf, + I dt I dxiicf, — 



(3.15) 



S ee S[u(x,t),u(x,t)} , (3.16) 
S<j, ee 5[u^(a;,t),u^(a;,t)] . (3-17) 



The additional term on the RHS of (|3 . 1 5|1 takes account 
of the non-inertial force due to the acceleration cb of the 
reference frame R' . 

The jacobian associated with the above transforma- 
tions is unity, as can be verified from the matrix elements 



0(xi,x 2 \ti,t 2 ) 



Su(xi,ti) 6u(xi,ti) 



Su <j) (x 2l t 2 ) Su < j > (x 2 ,t 2 ) 



5 ( xi - x 2 + / dt'<j)(t') ) 5(h - t 2 ) . (3.18) 

-T 



The operator which has the matrix elements (|3. 18[) can 
be written in any reasonable functional space of space- 
time dependent functions as 

6 = exp(J dt'(j)(t')-^\ . (3.19) 

The eigenstates of O are the momentum wavefunctions 
exp(zpx). Using a parity-preserving discretization of the 
Fourier space, the jacobian turns out to be 

det[0] = ]Jexp (ip J dt'<p(t') \ = 1 . (3.20) 



The Faddeev-Popov determinant A[u(0, t)] is also 
unity. In fact, consider the velocity field which has been 
"gauge fixed", with the help of a generalized galilean 
transformation, to u(0, t) — 0. We have, then, to sub- 
stitute the functional Taylor expansion of u^it) up to 
first order in <fr(t) in ()3.9j) . Defining g(t) = d x u(x, t)| x =o, 
we get 



«*(*)= / d£^)-(f>{t) + 0[<l> 2 {t)] (3.21) 

-T 



so that 



[u(0,t)] ee / D<j>8 



g(t) I dt'<P(t') - 4>{t) 

T 



= \det[S(t-t')-e(t-t')g(t)}\- 1 
Using, now, the identity 

det[X] = exp[Tr(lnX)] . 
we find, with A(t, t') = @(t - t')g(t), 

\n+l 



A[u(0,t)] = cxp 



V i '- Tr(A n ) 



71=1 



(3.22) 
(3.23) 

(3.24) 



since 



Ti(A n ) = J dhdt 2 ... dt n g(h)g(t 2 ) ... g(t n ) 
x0(ti - t 2 )Q(t 2 - t 3 ) ... G(t n -h)=0. (3.25) 

Above, we have used 0(0) = 0, which is the right pre- 
scription for the Heaviside function when the underlying 
stochastic differential equation is defined in terms of an 
Ito discretized time evolution [25| . Collecting all of the 
above results, we rewrite (|3.12[) as 



Dtp / Du^Du^u^t)} 



( [° d<p [ 

x exp I iS<f, + i J dt-r^ J dxticf, 



(3.26) 
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However, since u$(x, t) and u^(x, t) are integration fields, 
(j3.26[) becomes 



Z = J D(j> J DuDu6[u{0,t)] 
x exp (iS + i J dt^j- J dxu^j , (3.27) 
or, equivalently, integrating over 4>{t), 
Z= [ DuDu5[[ dxu(x,t)]5[u(0,t)]exp(iS) . (3.28) 



In view of (|3.27[) . we will assume, in all our subsequent 
considerations, that u(0,t) — 0. In other words, we have 
moved for good to the locally comoving reference frame. 
As fluid elements happen to stick (and spend a finite 
fraction of their times) in shock discontinuities, the lat- 
ter will be frequently hosted at the origin of the locally 
comoving reference frame. 

We stress, at this point, that the hypothesis of strong 
galilean invariance 19] is a fundamental ingredient here, 
since no role is given to the velocity boundary conditions 
at infinity, in the lagrangian formulation put forward in 
Eq. (13TT1) . 



IV. SHOCKS AND INTERMITTENCY 

The response functional (|3.28[) can be decomposed as 



Z = Z S 



(4.1) 



where Z s and Z s refer, respectively, to the cases where 
shocks and smooth field configurations are found at x = 
0, t = 0. Recalling the discussion of section II, it is clear 
that the statistical properties of large negative velocity- 
differences are all encoded in Z s . 

According to the probabilistic interpretation of Z, we 
note that Z s is not normalized to unit. Instead, Z s is 
normalized to the "intermittency factor" 7, where < 
7 < 1 is the fraction of time a shock is found at the 
origin of the locally comoving reference frame. 

Let u s {x) be a shock configuration, with velocity dis- 
continuity at x — 0, at its instant of creation (26j . Assum- 
ing that shock creation is uniformly distributed in time, 
one may write, for the probability density functional to 
get configuration uq(x) at time t = (see appendix A), 



x P[u 8 (x)]W[r],T;u Q (x),u s (x)] , 



(4.2) 



where Du s (x)P[u e (x)] is the probability measure for 
the creation of the shock u s (x), conditioned to be in 
the sample space of all shock creation events, and 
W[r], T; uq(x), u s (x)] is a weighting functional. 

Eq. (14. 2p is formally rigorous, but it is of hard 
practical implementation, due to the difficulty in 



getting information on the functionals P[u s (x)] and 
W[rj,T;uo(x),u s (x)]. Phenomenological arguments, 
however, can be helpful in order to replace ()4.2|) by more 
tractable expressions. 

Shocks are expected to have (i) mean interdistances 
of the order of the integral scale L and (ii) lifetimes of 
the order of T = L/U, where U is an estimate of the 
shock velocity jump. The prototypical Burgers shock is 
the stationary configuration 



(4.3) 



Even though (|4.3|) is a solution of the forceless Burgers 
equation, it can be used as a local approximation to gen- 
eral viscous shocks around the position of velocity discon- 
tinuity. Suppose that at time — T a configuration similar 
to (|4.3[) is created, and let g(U) be the probability den- 
sity that it has amplitude U. Due to property (ii) above, 
this shock is not going to be observed at time t = if 
T > L/U. The contribution to (|4.2[) provived by shocks 
with the local profile (|4.3|) is, then, estimated as 

r°° u r L / u r 

Z s = 7 J dUg{U)— J dTN J DiiDu 

x S[ / dxu(x,t)]5[u(0,t)]exp(iS) . (4.4) 



The velocity field u(x,t) in (|4.4[) satisfies the boundary 
condition 



u(x, — T) = u s (x; U) 



(4.5) 



We will work with the phenomenologically simplified 
result (14. 4|) in place of (|4.2p . In doing so, we conjec- 
ture that the asymptotic statistical properties of nega- 
tive velocity-differences are not affected by more detailed 
choices of shock parametrization. 

An interesting way to address the computation of (|4.4|) 
is to perform the corresponding path-integration over an 
appropriate subset of the functional space, which would 
consist of dominating configurations. That's precisely 
the purpose of the saddle-point method, applied in Refs. 
flU UjJ to the MSR turbulence context as a way to cope 
with the intermittency phenomenon. 

Saddle-point configurations, dubbed instantons, are 
associated with stationary values of the action. Taking 
functional derivatives of the MSR action (|3.4|) with re- 
spect to the integration fields, we get 



dtu + ud x u — vd^.u = —i / dx'u(x', t)D(\x — x'\) 



and 



dtu + ud x u + vd\u = 



(4.6) 



(4.7) 
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Also, when solving (|4.6p and (|4.7|) , we have to take into 
account the constraints [27J 



u(0,t) = , 
dxu(x, t) = 



(4.8) 
(4.9) 



Instanton solutions of (|4.6[) and (14. 7p . which hold for 
-T < £ < and satisfy l|£5j) , (gH) and (g^J) can be 
readily obtained: 



u(x, t) = u s (x; [/) , 
u(x, t) = . 



(4.10) 



It is worth mentioning that the solution for u(x,t) in 
(|4.10|) identifies Burgers shocks to instantons. Also, it is 
not difficult to find that the MSR action vanishes when 
evaluated for the fields given in (|4.10[) . As it is the stan- 
dard procedure in the saddle-point method, we expand 
the MSR action in a functional Taylor series around the 
instantons, retaining only quadratic fluctuations. We re- 
place, as a result, (|3.4p by 

S* = J dt J dxu[d t u + d x (u s ii) — vd^u] 

+— J dt J dxdx'u(x,t)u(x' ,t)D(\x — x'\) , 

(4.11) 

where the velocity boundary condition becomes, now, 
u{x,-T) = 0. 

In order to compute the pdf of negative velocity- 
differences, we introduce the characteristic function 

r°° tj r L /u r 

Z S {X) =7 / dUg(U)— I dTN / DuDu 



x6[J dxu(x, i)]S[u(0, t)} exp(iS* - iXz) , (4.12) 

where z is the velocity-difference evaluated at t = 0, 

z = -2U + u(C, 0) - u(-C, 0) . (4.13) 

The characteristic function Z S (X) can be exactly com- 
puted, in principle, since it is given in (|4. 1 2[) by a 
quadratic field theory. To evaluate Z S (X), the saddle- 
point method can be applied once again, this time in an 
exact way. The further saddle-points equations for u(x, t) 
and u(x, t) are 



dtu + d x (uu s ) — vd\u = 
= -i dx'u(x' ,t)D(\x - x'\ 



dtii + u s d x ii + vd\u = 

= X[S(x + C)-S(x-C)}S(t) 

supplemented by |gj5j| and (|43| . 



(4.14) 
(4.15) 



Observe that the viscosity term has the "wrong" sign 
in equation (|4.15p . To avoid the unbounded growing of 
u(x, t) for t > 0, we impose, as prescribed in Refs. [111 . 
[l6[ , the boundary condition u(x,0 + ) = 0. Integrating 
(I4.15P over the time interval [— e, e], with e — > 0, we get 
the "final condition" 



i(x,0-) = A[<y(ar-C)-*(a! + C)] • 



(4.16) 



Furthermore, we have the exact saddle-point result (see 
appendix B) 

S* - Xz = 2XU 

+ — J dt J dxdx 1 u(x , t)u(x' ,t)D(\x — x'\) . 

(4.17) 



It is interesting to note, due to (|4.17[) . that we do not 
have to worry in finding the specific solution for u(x, t). 
Eq. (|4.15p is solved, in the vanishing viscosity limit, by 



u(x, t) = X[6{x - x(t)) - 5{x + x(t))] , 



(4.18) 



where x(t) =C,-Ut Substituting fTTH]) into (j4T7]) . and 
taking £/£ << 1, we find 



S* - Xz = 
= 2XU + i 



DqL >2 
2U 



2UT/L 



dt[l 



(4.19) 



We get, from (|4~T2|) 

Z S (X) = 7 
x exp ^2iXU 
where 



dUg(U) / d V 
() Jo 

D Lc(ri) 2 



c(i?) 



dt[l 



(4.20) 



(4.21) 



-2 V 



The negative velocity-difference pdf is computed from 
the Fourier transform of the characteristic function, as 



1 

2^ 



dX exp(iXz)Z s (X) 



U 



T / iu ^ u) L d ^2^m 



x exp 



U(z + 2U) 



- 1 



2c( V )D Q L 



(4.22) 



Expression (|4.22|) gives, for large negative z, the asymp- 
totic pdf, 



(4.23) 
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where the dots refer to subleading contributions, and 

a = jLD g(0) f dr,c(r,) ~ Q.36jLD g(0) . (4.24) 
Jo 

The expression for the coefficient (|4.24l) is a testable pre- 
diction of the present theory. Alternative force-force cor- 
relation functions can be used to recompute (|4.24[) and 
compare it with the value to be found in further numeri- 
cal simulations. It is clear that in the eulerian framework, 
the intermittency factor 7 has to be replaced by 

7' = 2n£ , (4.25) 

where n is the number of shocks per unit length. 



V. CONCLUSIONS 

We have obtained, with the help of instanton tech- 
niques, the asymptotic form of the pdf of large neg- 
ative velocity differences in Burgers turbulence. The 
lagrangian picture of the Burgers flow was adapted to 
the MSR field theoretical framework, a procedure which 
proved to be an important technical improvement over 
the eulerian description. Lagrangian methods are, as a 
rule, welcome in the study of small-scale intermittency, 
since they cope in a natural way with the sweeping pro- 
duced by large scale motions. In the case of Burgers flow, 
sweeping produces shock advection, making it difficult to 
find out the statistical properties of velocity difference 
fluctuations. 

The introduction of lagrangian coordinates was carried 
out under the hypothesis of strong galilean invariance. 
We have found that the left tail pdf has the asymptotic 
form p{z) = a/|z| 3 , which agrees with the conjecture put 
forward in Ref. that this is so when strong galilean 
invariance holds. We have obtained an explicit expression 
for the critical amplitude a, which motivates the study 
of further numerical simulations of Burgers turbulence. 

Arbitrary shocks of the Burgers forceless model are 
identified to instantons, and taken, in the path-integral 
formulation of the response functional, as the dominant 
configurations for the determination of the velocity dif- 
ference fluctuations. We have bypassed the detailed clas- 
sification of all of these Burgers shocks at their creation 
events, by noting that relevant parameters of newborn 
shocks are their velocity jump, U, and extension, as- 
sumed to scale with the integral length L. Shocks are 
expected to have lifetimes of the order of L/U. In view 
of the role of the dimensional parameters U and L, we 
should regard the stationary Burgers shock (|4.10[) more 
as an illustration than as an essential ingredient in the 
formalism. 

We emphasize that the instanton distribution g(U) is 
not able to yield the left tail pdf of velocity differences on 
its own. Furthermore, the instanton contribution to the 
MSR action vanishes. The point is that the fluctuating 



background couples with the shocks and by the usual in- 
stability mechanism discussed in section II, large negative 
velocity differences are generated in the flow. One may 
wonder if this process of intermittency generation is anal- 
ogously found in the interaction between the background 
and coherent structures in Navier-Stokes turbulence. 
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APPENDIX A: DERIVATION OF EQ. (JOJ) 

Let us focus our attention on a flow which evolves un- 
der a particular realization of the stochastic force f(x, t). 
We also suppose that a random ensemble of initial con- 
figurations is given in the remote past (t — > —00), so that 
at any instant of time the possible velocity configurations 
yield a statistically stationary ensemble. 

We consider, now, the creation, at time — T, of a shock 
u s (x), localized at the origin of the locally comoving 
reference frame, which survives until time t = 0. Let 
T[u s , /] be the maximum value of T. Since the probabil- 
ity for the creation of the shock u s (x) in a time interval 
dT is also proportional to dT, we may write 

Z s = 7 J Du s (x)P[u s (x)] 

I r T[n 3 J] 

X ^ fl/ dTP[u o ,u s ;0,-T[u s J]}) f , 

1 [ u s, J\ Jo 

(Al) 

where P[uq, u s ; 0, — T[u s , /]] is the probability distribu- 
tion associated with the transition from the shock con- 
figuration u s (x), which was created at time — T, to the 
final configuration (at time t = 0) uq(x). For the sake of 
clarity, we note that 

P[u , u s ; 0, -T[u s , /]] = S[u (x) - L(x; [u s ,T, /])] , 

(A2) 

where L(x; [u s ,T,f]) is the velocity configuration which 
evolves from the shock u s (x) after the time interval T. 
Equation (lAljl can be rewritten as (14. 2|) . In fact, 

Z s = 1 I Du s (x)P[u s (x)} ^ *H f 1 dT 
J Jo V Jo 

x {S(r)-T[u s ,f])P[u o ,u s ;0,-T]} f 

= ran r dT f Dus{x) 

Jo v Jo J 

x P[u s {x)]W[r),T;u (x),u s (x)} , (A3) 

where 

W[ij,T; u Q (x),u s (x)} 

= (5( V - T[u s , f])P[u , u s ; 0, —T])f . (A4) 
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Observe that if T[u s , f] does not depend on f(x,t) 
(which may be a useful approximation), then 

W[r],T;u (x),u s (x)} 

= 5(r)-T[u s ])(P[u ,u s] 0,-T]) f 

(A5) 

and, therefore, 

r i fT[u s ] 

Z s = 7 / Du s {x)P[u s {x)]— - / dT 
J T[u s \ Jo 

x (P[u ,u s ;0,-T}) f . (A6) 

Above, the averaged probability (P[u , u s ; 0, — T])/ can 
be given as the MSR path-integral expression (|3.3[) . 

APPENDIX B: DERIVATION OF EQ. (|4lTl) 

There is some subtlety in the saddle-point evaluation of 
characteristic functionals like (|4.12|) . Since we are con- 



sidering in (|4.12|) the evolution up to time t — 0, one 
could object that the u(x,0 + ) = boundary condition 
sounds too loose. Actually, in order to apply the saddle- 
point method to (|4.12[) . the time evolution is extended 
to t — > oo. Saddle-point solutions are, then, such that 
u(x,t — > oo) = u(x,t — > oo) = 0. Note that the time 
extension does not change the value of Z S (X), once ve- 
locity configurations are integrated out at t — > oo in the 
path-integral ()4.12[) . 

Taking these remarks into account, we multiply both 
sides of Eq. (|4.15|) by u(x, t) and integrate them over 
space and time. We find 

/ dx J dtu[dtu + u s d x u — vd^u] 
= A[«(C,0)-«(-C,0)], (Bl) 

where we have used the boundary conditions u(0, x) = 0, 
u(x, t > 0) = and u(x, —T) = 0. Eq. (|4"T?1) follows 
straightforwardly from the substitution of (|B1[) in the 
expression for S* — Xz, as it enters into (|4.12p . 
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